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Abstract 

An anomaly in time distribution of neutrinos from the ISIS pulsed beam 
stop source observed by the KARMEN collaboration is discussed. We show 
that the anomaly can be interpreted as a superposition of two exponentials, 
both having time constants consistent with the [i + lifetime of 2.2 [is. It 
is assumed that they both originate from muon decays at rest. One of 
them describes the time distribution of the prompt neutrino events, while 
the other describes the time distribution of events from delayed decays of 
slowly moving {(3 ~ 0.02) particles in the KARMEN calorimeter. 

We propose here that these particles are produced in exotic decays of 
positive muons u + — > e + + X, resulting in the second exponential time 
distribution shifted by the time of flight with respect to the time distribu- 
tion of neutrino events. This model gives an aceptable fit to the KARMEN 
data if X has a mass of 103.9 MeV. The possible decay modes of this new 
massive neutral particle are discussed. This hypothesis can be experimen- 
tally tested in the near future by studying the low energy part of the e + 
spectrum in the /i + decays. 

1 Introduction 

The KARMEN collaboration, which studies the interactions of neutrinos from 
the stopped n + decay chain at Rutherford Appelton Laboratory, has reported 
an anomaly in the time distribution of neutrino events in the KARMEN 
calorimeter. There is a bump of events in the distribution which was expected to 
be well described by a single exponential with a time constant equal to the muon 
decay lifetime. 

It was suggested that a new weakly - interacting half spin particle with a 
mass 33.9 MeV/c 2 from initial 7r + — > fi + + x decay , which travels with velocity 
/9 = f/c~l/60 and decays in the KARMEN calorimeter after a mean flight path 
of 17.5 metres, can be a possible candidate for an explanation of this anomaly. 
Since the visible energy in the calorimeter was in the range from 11 to 35 MeV 
it was also assumed that these new particles deposited only a fraction of their 
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full energy in the calorimeter. The hypothetical 7r + — > fi + + x decay has been 
searched for at PSI by studying muons from pion decay in flight 0. The result 
gives a branching ratio of less than 7 x 10~ 8 at 95% Confidence Level (C.L.) for 
this decay mode leaving only a narrow allowed region in the branching ratio if 
combined with the results of ref. ||. 

Another possible interpretation of the anomaly proposed by the KARMEN 
collaboration is related to the direct production of a new particle in the proton 
collision with the neutrino target. However, the collaboration has not yet stated 
whether it has devised one. 

Concerning the possible explanation of the anomaly by exotic muon decays 
it was concluded that "...the narrow clustering of the anomalous events rules 
out contributions from fi + decay, which spreads over a much longer time interval 
(compared to the 7r + - life time of ~ 20 ns) up to 10 /is after beam-on target." 

!• 

This conclusion is the motivation for present letter, where we show that the 
KARMEN anomaly could be explained by contributions from the exotic fi + de- 
cays. 

The paper is organised as follows. In section 2 we analyse the time distribution 
obtained by the KARMEN collaboration and give a possible explanation of the 
anomaly. As an example we consider the fi + — > e + + X decay into a positron and 
a new boson with the mass of 103.9 MeV/c 2 . In sections 3,4 the possible X-boson 
phenomenologies and experiments on searching for the rare decay /x + —> e + + X 
are discussed, respectively. Section 5 contains concluding remarks. 

2 Analysis of the neutrino time spectrum 

Figure 1 shows the time distribution found by KARMEN with respect to beam- 
on target events. Analysis of this time distribution was performed in ref. [0] after 
subtraction of the cosmic ray induced background, which was found to be constant 
and equal to 60 events per 0.5 

This spectrum was fitted in ref. [|TJ by a superposition of Gaussian and expo- 
nential signals from muon decay neutrinos having a 2.2 /is decay constant. The 
fit procedure performed with the minimization routines of MINUIT Q results 
in x 2 of 9.7 for 14 degrees of freedom (DF) corresponding to 79% C.L. for this 
hypothesis. Gaussian, which is centred at (3.6 ± 0.25)/is with a full width at half 
maximum(FWHM) of (1.3 ± 0.5)/xs consisted of 83±28 events associated with 
the signals from decay in flight of the proposed particle in the fiducial volume of 
the KARMEN calorimeter. 

We tried to reproduce these results by using the same MINUIT minimizing 

3 The experimental points were reproduced from the original KARMEN publication jij . The 
error bars were calculated taking into account the errors due to subtraction of the background. 
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routines. Our results and results of the KARMEN fit were found to be in a good 
agreement: our fit procedure yielded \ 2 °f 9-8 f° r 14 degrees of freedom (DF) 
corresponding to 77% CL. for this hypothesis. Gaussian was found to be centred 
at (3.7 ± 0.26)/xs with a FWHM of (1.3 ± 0.5)/is and with 85±32 events. The 
result of our fit procedure is shown in Figure 1. 

Having this good agreement in mind, we then performed an exponential fit 
of the KARMEN time distribution assuming that this spectrum is described by 
a sum of two exponentials: the first exponential describes the time distribution 
in the interval 1 - 3.3 /is , and the second one describes the time distribution 
in the interval 3.3 - 10 fis. Figure 2 shows the result of the fit (solid lines). The 
MINUIT minimizing routines yield a x 2 of 9.7 for 15 DF corresponding to 84% 
C.L., and time constants of (2.29 ± 0.34)/xs and (2.1 ± 0.6)/xs for the first and 
the second exponentials respectively. Note that there is no significant difference 
in Confidence Levels of our and KARMEN 's hypotheses. 

The time width of the transition region between two exponentials is rather 
small and corresponds to ~ 0.5/is. We try to fit this region more precisely using 
a Gaussian shape of the transition curve instead of a step-function. However, no 
significant difference between two results has been found. 

A similar description of the time distribution was obtained by superimposing 
an exponential describing the distribution of events in the time interval 3.3 - 
10 jis upon an exponential signal having a 2.2 jis decay constant for the full 
time interval. Figure 3 shows the corresponding time distribution of events after 
subtraction of the exponential part with 2.2 jis. A direct exponential fit of the 
time distribution of events in this region gives a time constant of 2.26 jis close 
to the value of muon decay constant; however the error is rather large because of 
poor statistics. 

Since the time constants were found to be equal within the error limits to 
the muon decay constant we may assume that all events in the KARMEN time 
distribution originate from the muon decays. 

Thus, our interpretation of the KARMEN anomaly is the following. Positive 
muons from 7i + decays stop in the target ( muon range in the target is ~ 1 cm) 
and decay giving a flux of neutrinos which decreases in time exponentially with 
the muon decay constant. The first exponential shown in the time spectrum in 
Figure 2 is related to the interactions of this flux of neutrinos in the KARMEN 
calorimeter. An excess of events after ~ 3.6 jis shown in Figure 3 is related to the 
decays of slow moving particles produced in muon decays. Therefore, the flux of 
these particles and fraction of their decays in the fiducial volume of the KARMEN 
calorimeter also decrease in time exponentially with the muon decay constant, 
resulting in the second time distribution but shifted by the time of flight with 
respect to the time distribution of neutrino events. Since the time width of the 
transition region between two exponentials is rather narrow (see Figure 2), we 
suggest that these new particles originate from the two-body decays of positive 
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muons, namely /i — > e + + X. We note that the RMS of fluctuations of time needed 
for X bosons to pass from the neutrino target to the calorimeter and decay inside 
its volume was estimated to be of the order of 0.2 /j,s. The effect of the convolution 
of these fluctuations and the exponential from muon decay on the exponential 
shape of this curve and its decay constant was found to be small. 

The mass of the X boson calculated for (3 X = 1/60 is found to be (103.9±0.1) 
MeV/ c 2 , where the error corresponds to the spread of X's velocity estimated from 
the time width of the transition region between two exponentials. The excess of 
events after subtraction of the contribution from the neutrino was found to be 
AN = 108 ± 41 events. 

The limit on branching ratio Br = — > e + + X)/T(fi — > all) as a function 
of the X-boson lifetime can be recalculated from the KARMEN result for n + — > 
fi + + x decay mode taking into account the difference in the number of excess 
events found after the fit procedure. The result is shown in Figure 4. 



3 Possible X-boson phenomenologies 

Suppose that X-boson is a neutral scalar particle. The Lagrangian describing the 
fi — > e + X decay has the form 

L = X[h lxeL jl L e R + h IMeR fi R e L + /i.e.] (1) 

The /i — > e + X decay width is given by the formula 

r( ^ e+ X)^(^ + ^)(l-^). [( l + ^f-4|]^ (2) 

For a non-relativistic X-boson with a speed vx = ^jc we find that its mass is 
rax = 103.9 MeV and the fi — > e + X branching ratio is 

BrQM - eX) » 0.4 • 10 13 [^ ei + h% R } (3) 

For the interaction (1) the X-boson decays via the virtual \x meson X — > (//* — > 
ez/z/) + e and its decay width is 

T(X -» e+e^P) ~ 0((7r 2 r(^ -» e^)(^ eL + ^)) (4) 

As it follows from the observed number of the anomalous events due to the 
X-boson decay in the KARMEN experiment, the ratio 

K = Br(/j, — > e + X)T(X — ► vis) (5) 



is 



K ~ 0(3 • lO- 11 ^- 1 ) (6) 
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From the formulae (2 -6) we find that Br(fi — > e + X) ~ O(10~ 2 ) and r(X — > 
\ievv) ~ O(10 9 )s, that is probably at the level of the contradiction with the 
experimental bound Br(/i — > e + X) < 3 ■ 10~ 4 which is valid however for 
mx > 103.5 MeV (see Section 4). It is possible to invent a model where the 
Br(p, — > e + X) is less than O(10 -4 ). Let us introduce the hypothetical heavy- 
charged lepton L and consider the decay of the X-boson through the virtual 
charged L- lepton 

X -> (L* -» ei/P) + e (7) 

For the L-lepton mass mi = 100 GeV, the coupling constant of L- lepton with 
electron and X-boson /iL e x = 0.1, the L-lepton mixing with electron e 2 eL = 10~ 2 , 
we find that Br(fi -> eX) = 0(1O~ 5 ). So, if our interpretation of the KARMEN 
anomaly is correct we expect a rather large branching ratio Br(fi — > e + X) > 
O(10 -5 ). Therefore, it would be very interesting to perform an experiment on the 
search for // — > e + X decay with improved sensitivity. These scenarios predict 
the X-boson lifetime t(X) > O(10 6 ) s that probably leads to problems with 
astrophysical and cosmological constraints. Moreover, these scenarios predict an 
anomalously high X-boson decay energy < E V i s >~ 50 MeV compared to the 
observed visible energy Eyis ^ 

(11 - 35) MeV. 

To avoid these problems, let us consider the model with two additional neutral 
scalar bosons <pi and <p2 and suppose that the 0i-boson is lighter than the X-boson 
and the 02-boson is heavier than the X-boson. The Lagrangian describing the 
X-boson decay has the form 

L = XX(pi(f)2 + h e <p 2 ee + h v $\i>v. (8) 

The X-boson decays via the reaction 

X ^ + (0* _> e+e") (9) 

through the virtual 02-boson and the neutral </>i-boson decays into neutrino an- 
tineutrino pair. For such a scenario it is possible to have < T vis >~ 30 MeV and 
r(X) < 10~ 2 s. For instance, for m x = 103.9 MeV, m 01 = 73 MeV, m 02 = 1 
GeV, h e = 10 -3 , A = 10~ 2 m^ and hp = 10~ 6 we find that the average visible en- 
ergy in the X-decay is < T ms 26 MeV and r(X) ~ 3 ■ 10~ 3 s, r(0 x ) ~ 2 • lO" 10 
s, Br([i —> e + X) « 10 -13 , r(02) ~ 3 • 10 -17 s. So in this rather artificial model 
all particles are not long-lived and it is possible to have < T vis > in agreement 
with KARMEN data. 

Note also that it is possible to have a model where the 02-boson is the Higgs 
boson. For instance, for the Higgs boson mass = 100 GeV and A = we 
find that r(X 0ie+e~) ~ 0.7 ■ 10 7 s and Br(/2 -»■ eX) ~ 2 ■ 10^ 4 . In this 
case to have r(X) < 10~ 2 s we have to postulate an additional interaction of the 
X-boson with neutrinos. 
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4 Experimental search for the fi + —> e + +X decay 

Since X's are weakly-interacting and relatively long lived particles, their flux 
would penetrate any type of calorimeter without significant attenuation and can- 
not be observed effectively in the detector via its decays. A more effective exper- 
imental signature for fi + — > e + + X decays would be a peak from mono-energetic 
positrons in the low energy part of the positron spectrum corresponding to the 
principal mode of /i + decays. 

Direct experimental searches for a peak from fi + —> e + + X decays were 
performed in a few experiments f|. For the lowest positron energy region from 
1.6 to 6.8 MeV, corresponding to the mass region of 103.5 - 98.3 MeV/c 2 , the 
branching ratio limit of — > e + + X)/T(fi + — » all) < 3 x 10~ 4 was obtained 
in ref. § from the analysis of published data on muon decays in the hydrogen 
bubble chamber obtained by Derenzo 0. The result is limited by a poor ( ~ 
30-40%) momentum resolution in this region. 

It would be interesting to search for a peak in the low energy part of the 
positron spectrum by using, for example, a high resolution Ge-detector ( Nal 
detector) with a typical energy resolution of the order of KeV (~ 100 KeV). 
The low energy muons could be stopped in the active volume of the detector. In 
this case the low energy spectrum of the positrons could be measured without 
significant distortions. 

The estimate of the rate from the principal muon decays shows that, using 
Ge-detector, one could expect a sensitivity to the fi + —> e + + X decay branching 
ratio of better than 10~ 6 . Here it is assumed that the intensity of stops is of the 
order of 10 4 /i + /s and that the exposure to the muon beam is ~ 1 month. 

5 Conclusion 

The anomalous time distribution of the events from the KARMEN detector can 
be acceptably described by a superposition of the time distribution of neutrino 
signals following the expected 2.2 /j,s muon decay constant and a weaker exponen- 
tial signal from delayed decays of the new massive neutral particles also originated 
from the muon decay at rest and shifted in time due to their time of flight. 

It is proposed that a new boson X with the mass of 103.9 MeV/c 2 which 
might be produced in the decays of positive muons /z + — > e + + X could be a 
possible candidate for an explanation of the anomaly. The visible energy from 
the possible decay modes of X's can be within the allowed region from 11 to 35 
MeV/c 2 . This hypothesis can be experimentally tested in the near future at a 
level of sensitivity of better than 10~ 6 . 

It should be noted that the original interpretation of the KARMEN anomaly 
as a new neutral fermion x produced in 7r + — > [i + x decay with a mass m x = 33.9 
MeV has been discussed in ref. 0. It has been shown that a mainly sterile 
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neutrino scenario is compatible with all laboratory constraints, although there 
are problems with astrophysical and cosmological constraints. Our interpretation 
with the scalar X-boson produced in the fi + — > e + X decay is in the spirit of the 
familon scenario ||. 
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Figure 1: Time distribution of events in the KARMEN calorimeter after the 
subtraction of the cosmic background? The data are fitted to an exponential with 
the 2.2 lis decay constant on which is superimposed a Gaussian signal centred at 
3.7 jis. The fit procedure results in x 2 of 9.8 for 14 degrees of freedom. 
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Figure 2: Time distribution of events in the KARMEN calorimeter after the 
subtraction of the cosmic background? The solid curves are a fit to the points by 
a sum of two exponentials. The first exponential describes the time distribution in 
the region from 1.0 to 3.3 us and the other in the region from 3.3 to 10 fis with 
time constants of (2.29 ± 0.34)yus and (2.1 ± 0.6)fis, respectively. The broken line 
corresponds to the extrapolation of the first exponential. The fit procedure results 
in x 2 of 9.7 for 15 degrees of freedom. 
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Figure 3: Time distribution of excess events left after subtraction of the exponen- 
tial component with 2.2 fis decay constant. The solid line is the exponential fit to 
the points in the interval 3.3 - 10 /is which gives a time constant of 2.26\xs. 
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Figure 4: The branching ratio — >■ e + + X) /T(fi + — > all) as a function of the 
lifetime of the X -boson calculated from the number of observed events ( AN = 91 
events) using the analogous KARMEN plot for a decay mode of n + — > ji + + x . 
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